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Abstract 

Preliminary ab-initio applications of many-body Green's functions theory to 
the ground state of ^He suggest that high accuracy can be achieved in the 
so-called Faddeev-random-phase-approximation method. We stress the poten- 
G^ \ tialities of this approach for microscopic studies of medium-large nuclei and 

^ ■ report applications to IsOd and lpO/-shell nuclei. In particular, we discuss the 

psj . role of long-range correlations on spectroscopic factors and their dependence 

Q_il on asymmetry. 

D 



1 Introduction 

The advent of radioactive ion beams has opened the possibility of reaching unexplored regions of the 
nuclear landscape. This has lead to the realization that single particle levels evolve with changing the 
number of neutrons or protons [1,2], sometimes invalidating the conventional magic numbers known for 
stable isotopes. Another unexpected observation is that spectroscopic factors for the same states (which 
we will also refer to as quasiparticles) may change dramatically with proton/neutron asymmetry [3]. 
^ ■ This implies non trivial (and not yet well understood) changes in nuclear correlations at the driplines. 

This talk reports about on going investigations of the evolution of quasiparticle properties using large 
^ \ scale microscopic calculations. 

^ \ From the point of view of ab-initio nuclear structure, great progresses have been made for light 

^ ' nuclei. There exist now a wide range of techniques that solve the few -body problem exactly [4]. Moving 

■ to medium and large nuclei is made difficult by the exponential grow in degrees of freedom. In this 

^ ! case, techniques like coupled-cluster (CC) [5-7] and self-consistent Green's functions (SCGF) [8-10] 

are preferable. For both methods, the number of configurations to be dealt with scales more gently 
^ ' with the increasing particle numbers and the infinite summation of connected diagrams guarantees size 

O ■ extensivity and accuracy. Applications to medium-heavy nuclei are (for the moment) limited to systems 

around closed subshell structures. 



^ \ For open shell systems, standard shell model calculations remain the best option. However, one 

■ faces questions about what effective interactions and charges have to be used when extrapolating to the 

driplines. One could attempt to bridge shell model Hamiltonians to realistic nuclear forces by employing 
the above ab-initio approaches and studying quasiparticle states and their interactions for different shell 
closures across the nuclear chart. In this respect, the SCGF is the method of choice since the one-body 
Green's function already contains, in its energy representation, all information of quasiparticle properties. 
As we will note below, self-consistency (SC) helps in improving accuracy and implies the fulfillment of 
conservation laws. Moreover, the strong connection between the Green's functions formalism and the 
experimental data allows to gain unmatched insight into the nuclear dynamics. 

The state of the art SCGF calculations employ the so-called Faddeev random phase approximation 
(FRPA) to expand the nuclear self-energy. This method is detailed in [9, 11, 12] and is summarized di- 
agrammatically in Fig.[T] The Ii{uj) and g^\u>) propagators carry information on the collective motion 
of particle-hole (giant resonance) and two-particle or two-hole (pairing like) configurations. The interac- 
tion between nucleons and collective excitations of the nucleus [11, 12] are accounted for by resumming 
infinite series of Feynman diagrams, like the one depicted in Fig. [it. The inclusion of RPA phonons 
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Fig. 1: The self-energy separates exactly into a mean field term and the polarization propagators R{u!) for 

the 2plh/2hlp motion, as shown in a). The double lines represent fully correlated propagators. Upon expansion 
of R{ui) in Feynman diagrams, one obtains the series of diagrams b) for the self-energy. The diagram c) gives an 
example of the contributions to R'^p^'''{uj) that are summed to all orders by the FRPA method. 



makes this method particularly suited for medium-heavy isotopes, where collective states are important. 
Nevertheless, we will show in Sec. |2] that rather accurate results are possible even for nuclei as small 
as '^He. We will further discuss the accuracy and convergence of quasiparticle energies in Sec. [3l And 
finally review ongoing applications to spectroscopic factors in IsOd and IpOf shell systems, including 
cases of large proton-neutron asymmetry. 

2 Ab-initio Calculations of ^He 

The spectral function of '^He has been computed in an harmonic oscillator basis using the Viowk in- 
teraction [13, 14] (derived from the Argonne vl8 potential with a sharp cutoff of A= 1.9 fm^^). We 
solved for the FRPA self-energy including all contributions to third oder, as explained in Ref. [12], and 
used the Koltun sum rule to extract the binding energy [10, 15]. The calculations employed the intrinsic 
Hamiltonian, Hint = T + V — T^m, to remove the contribution of the kinetic energy of the center of 
mass. 

Figure |2] shows the convergence of the total energy with respect to the oscillator frequency h^l and 
the basis size. Up to 17 shells iN,max=^6) and partial waves with l< 5 were used. For the largest space, 
results are basically converged (to within 10 keV) and approach exponentially to -29.00 MeV. This has 
to be compared with the exact Faddeev-Yakubowsky result of -29.19(5) keV [5]. An estimate of the 
self-consistency effects — given by the dashed line in the right panel — further lowers the total energy and 
brings it closer to the exact result. At the time of writing, the converged result for full self-consistency is 
not known and more work is in progress to answer this question [16]. 

We note that, based on the perturbative analysis of Ref. [17], the third order algebraic diagram- 
matic construction method [ADC-(3)], the similar «o?i-self-consistent FRPA, and coupled-cluster theory 
with singles and doubles (CCSD) are all expected to give similar results for the binding energy. Do- 
ing self-consistency in Green's functions theory introduces implicitly several other diagrams. Among 
these are the same fourth and fifth order diagrams evaluated by the triple corrections to the CC method 
[CCSD(T)]. This is indeed the picture that can be inferred from comparing Fig. |2] and the analogous CC 
study of Ref. [5]. Further confirmation comes from analogous studies of atomic systems [18]. 

Self-consistency introduces effects of correlations that are beyond the bare FRPA level and are 
expected to improve sensibly the accuracy of results. Following this approach has two main advantages. 
First, the corrections are included in the self-energy and therefore apply to all the quantities derived from 
it. These include the total energy, quasiparticle energies, overlap wave functions, scattering observables 
and so on. Second, it is known from the theorems of Baym and Kadanoff [19, 20] that this procedure 
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Fig. 2: Convergence of the FRPA binding energy for ^He. Left: dependence on the oscillator parameter and 
basis size (without self-consistency). Right: convergence for fixed hil=l4 MeV (black line) and corresponding 
correction due to approximate self-consistency (dashed green line). The horizontal lines mark the exact Faddeev- 
Yakubowsky result. 



selects those additional diagrams that allow to satisfy general conservation laws. 



3 Calculations of Quasiparticle States 

The self-consistent FRPA approach has been employed in Ref. [8] to study the single particle spectrum 
around the Fermi surface of ^^O. A large oscillator basis of eight shells was employed. The resulting 
energies for the addition or removal of a nucleon are shown in Fig. [3] as a function of the oscillator 
length. The two panels show the results obtained for the regularized Vucom [21] interaction and the 
Argonne vl8 [22] potential (both without the Coulomb force). For the sole case of the Argonne potential, 
the G-matrix technique was employed to resum ladder diagrams from outside the model space. This 
approach allows to properly include the effects of short-range correlations [9,23]. While quasiparticle 
and quasihole energies obtained with Vucom still carry some dependence on the oscillator parameter, the 
corresponding spin-orbit slitting appears to be rather well converged. For the Argonne vl8 calculation, 
the dependence on the oscillator parameter is sensibly weaker, due to the inclusion of ladder diagrams 
from outside the model space. 

The results for spin-orbit splittings of the Op and Od orbits are reported in Tab. [T] and compared 
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Fig. 3: Neutron quasiparticle energies for ^^O, as obtained for the Vucom (left) and the Argonne vl8 (right) 
interactions [8]. For the latter, a G-matrix resummation of short-range diagrams was employed. 
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Fig. 4: Left: neutron quasiparticle energies for ^^Ni obtained for the N3LO interaction as a function the the oscilla- 
tor frequency. A monopole correction was added to constrain the ph gap to its experimental value at hU=lQ MeV. 
Right: correlation between the spectroscopic factors for the lp3/2 and O/7/2 orbits and the ph gap. Each line 
correspond to different HQ and sizes of the model space and is obtained by varying the monopole correction term. 



to the experimental values for neutron single particle states (as observed in O and O). The UCOM 
renormalization procedure tames the short-range correlations at the level of two body forces and en- 
hances the non-locality of the interaction. This leads to larger spin-orbit splittings than the vl8 case and 
closer to the experiment. However, the energy gaps between quasiparticle energies and the Opi/2 hole 
are overestimated by the Vucom interaction. This is a consequence of the well know issue of nuclear 
radii, which are usually underestimated by this interaction without three-nucleon forces [25]. The total 
energy of ^^O with this interaction has also been studied more recently, using similar model spaces [26]. 
Note that the Op splitting of 3.1 MeV obtained for the Argonne vl8 potential is very close to the vari- 
ational Monte Carlo result of 3.4 MeV (which was derived for the older vl4 model) [24]. This gives 
further confirmation of the accuracy of our calculations. The remaining discrepancy with respect to the 
experiment should not be attributed to the many-body method but to the neglect of three-nucleon forces 
in the Hamiltonian. The vl8 interaction with a G-matrix treatment of the short-range physics reproduces 
particle-hole (ph) gaps in much better agreement with the observation. 

Recent improvements in the numerical codes for calculating the FRPA self-energy have allowed 
to extend the above calculations to larger nuclei in the Oplf shell. For ^^Ni, calculations were performed 
in spaces up to ten oscillator shells [9]. In this work, the chiral N3L0 interaction [27] was modified 
by a simple correction of the monopole interaction. This correction is needed with most realistic two- 
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Table 1: Spin-orbit splittings and ph gaps of ^^O (in MeV) obtained for the Vucom and Argonne vl8 interac- 
tions [8]. A G-matrix resummation of ladders outside the model space is employed for vl8. The experimental 
values refer to the spectra of the corresponding states in ^^O and ^^O. 
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Table 2: Energies (in MeV) and spectroscopic factors (as a fraction of the independent-particle model) for tran- 
sitions to the IpO/ valence orbits of ^^Ni, obtained for Ml — 10 MeV [9, 16]. The third and fourth columns 
correspond to the sole contribution of correlations from outside the model space (SRC) and to the full FRPA result. 
The corrections AZ^ are obtained by comparing shell model and FRPA in the pf shell alone and are added in the 
sixth column. The experimental single-particle energies are taken from [30]. The measured spectroscopic factor 
for transfer between the ground states of ^^Ni and ^^Ni is from Ref. [35]. 

nucleon Hamiltonians in order to reproduce correctly the experimental gap between the 1^3/2 and the 
0/7/2 orbits [28]. The dependence of single particle energies on the oscillator frequency, Ml, is given 
by the left panel of Fig. |4l The results show limited deviations for a large range of frequencies. Part of 
the residual dependence is also due to the prescription used for the monopole correction, which is also 
dependent on Ml. Interestingly enough, it is found that the calculated spectroscopic factors for these 
single particle states are only sensible to the ph energy gap. This is shown by the right panel of Fig. |4j 
where calculation are reported for different sizes of the model spaces and monopole corrections [9]. It 
follows that, even if calculations are not perfectly converged, stable predictions for the spectroscopic 
factors can be obtained by constraining the ph gap. A study of the contributions to the spectroscopic 
factors from different types of correlations is reported in the next section. 

4 Results for Spectroscopic Factors 

Spectroscopic factors for particle and hole states around closed shell nuclei are of particular interest since 
they carry important information regarding nuclear structure [29]. Strong deviations from unity signal 
the onset of substantial correlation effects and imply the existence of non trivial many-body dynamics. 
Understanding how spectroscopic factors change when moving close to the drip lines is important to 
constrain theoretical models of radioactive isotopes. 

4.1 Role of Long-Range Correlations 

Spectroscopic factors for the quasiparticle and quasihole peaks around ^^Ni are reported in Tab. [2] [9, 16]. 
These correspond to the single particle levels displayed in Fig.|4]for Ml=\0 IVIeV and were calculated for 
the same monopole corrected N3L0 interaction (which reproduces the ph gap at Z,N=28). Theoretically, 
spectroscopic factors are defined as the norm of the overlap function V'a(r) for the nucleon transfer to a 
given quasiparticle state a. This can be obtained from the nuclear self-energy by the relation [10,29] 



where S^^(a;) =< V'q,|S(cj)|^/;q, > is the expectation value of the self-energy calculated for the overlap 
function itself but normalized to unity (/ iV'ap = 1). 

The third column in Tab. [2] shows the results obtained from the self-consistent self-energy cal- 
culated by including in Eq. ([Hi only the contributions induced by degrees of freedom from outside the 
model space. As the N3L0 interaction is rather soft compared to other realistic models the effects of 
short-range correlations (SRC) are small. For interactions such as the Bonn or Argonne models, one 
should expect a quenching up to about 10% or more [29,31]. These early results are confirmed by recent 
electron scattering experiments [32-34]. The complete FRPA results, in the model space of 10 oscillator 
shells plus SRC, are given in the fourth column. For the transition between the ^^Ni and ^^Ni ground 
states our result agrees, within the error bar, with the experimental value of Ref. [35]. The remaining 
orbits have similar quenching and are in line with values expected for stable closed-shell nuclei. There- 
fore, for this nucleus the dominant mechanisms in reducing the spectroscopic strength is to be looked for 
in long-range correlations, involving configurations among a several shells around the Fermi surface. 

We now turn to the question of what are the effects of correlations induced by configurations 
beyond 2plh and 2hlp. Since extra ph excitations are induced by the self-consistency approach one 
can expect the SC-FRPA to be accurate for closed shell nuclei, even though one does not perform a 
full configuration mixing. To show that this is the case, we performed both SC-FPRA and shell model 
calculations in the IpOf shell [36]. This was done using the same interaction and single particle orbits 
employed in the large oscillator basis. The ANTOINE shell model code was used [28, 37]. As an 
example, the spectroscopic factor for the ground state of ^^Ni, was predicted to be 0.82 by FRPA 
and 0.79 by the shell model with up to 6p6h states. Comparing with the full space result of 0.65 in 
Tab. 121 one can infer that about half of the quenching is driven by degrees of freedom from outside 
the pf space. The difference between shell model and FRPA, AZq,=0.79-0.82=-0.03, is rather small as 
expected. The fifth and sixth columns in the table show the corrections AZq, for all single particle orbits 
and the resulting spectroscopic factors when the AZqS are subtracted from the full space FRPA results. 

It remains clear that the extra correlations from the shell model (i.e., not included by FRPA) will 
become crucial in open shell nuclei. Nevertheless, these results give confidence that calculations based 
on Green's functions theory can lead to highly accurate results for quasiparticles in closed shells. This 
opens a new way to derive the properties of shell model Hamiltonians from realistic forces. 

4.2 Proton-Neutron Asymmetry 

Understanding how spectroscopic factors change when moving close to the drip lines is important to 
constrain theoretical models of radioactive isotopes. First information on these features has recently 
become available using one-nucleon knockout experiments in inverse kinematics. In general, it is found 
that spectroscopic factors do change with proton-neutron asymmetry and the quenching of quasiparticle 
orbits close to the Fermi surface become stronger with increasing separation energy [3]. 

Figure [5] shows FRPA results for the spectroscopic factors of quasiparticles around ^^'^Sq and 
AOfiOQ^ [38]. These are based on the realistic chiral N3L0 interaction and employ also a G-matrix 
to account for the effects of short-range conelations. A dependence on proton-neutron asymmetry is 
indeed observed in the FRPA, with the spectroscopic factors becoming smaller with increasing nucleon 
separation energy. A dispersive optical model analysis, which is constrained to data up to ^'^Ca, has also 
been extrapolated to neutron rich Ca isotopes with similar findings [39]. However, for both analyses the 
change in magnitude is significantly smaller than the one deduced from direct knockout data [3]. 

We stress that the calculations of Fig. [5] fail in reproducing energy gaps across shells and the 
corresponding location of giant resonances. This is a typical finding for soft low-energy interactions 
(see, for example, Vucom in Tab. [T]) which are incomplete without multi-nucleon forces. Based on 
the findings of Sec. 14. 1[ collective excitations from Fig. [T]: are the most important degrees of freedom 
governing the reduction of spectroscopic factors [16,29,38]. These are properly accounted for by the 
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Fig. 5: Spectroscopic factors obtained from partially self-consistent FRPA. All numbers are given as a fraction of 
the independent particle model value and refer to transitions from ground state to ground state. The points refer to 
knockout of a nucleon from the isotope indicated nearby. The lines are a guide to the eye. 



FRPA approach. However, either realistic three-nucleon forces or phenomenological corrections will be 
needed to achieve reliable predictions. It is plausible that the dependence on asymmetry seen in Fig. [5] 
may become more substantial once FRPA calculations with improved interactions will be available. 

5 Conclusions 

We reported about ongoing calculations of the binding energy of ^He using the FRPA method within the 
framework of self-consistent Green's functions theory. The results show that this approach can achieve 
high accuracy and provides a promising method for ab-initio calculations in medium-heavy nuclei. As 
an example of its application, we discussed quasiparticle states in ^^'^^O, ^'^'^'^Ca and ^^Ni and the N/Z 
dependence of the spectral strength. Spectroscopic factors of nuclei around ^^Ni are found to receive 
very small corrections from missing many-body correlations. This approach opens a new way to derive 
the properties of shell model Hamiltonians from realistic nuclear forces. 
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